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Abstract. We have conducted a deep survey (r.m.s noise ~ 17/iJy) with the Very Large Array (VLA) at 1.4 GHz, 
with a resolution of 6 arcsec, of a 1 deg^ region included in the VIRMOS VLT Deep Survey. In the same field we 
already have multiband photometry down to Jab ~ 25, and spectroscopic observations will be obtained during the 
VIRMOS VLT survey. The homogeneous sensitivity over the whole field has allowed to derive a complete sample of 
1054 radio sources {5a limit). We give a detailed description of the data reduction and of the analysis of the radio 
observations, with particular care to the effects of clean bias and bandwidth smearing, and of the methods used 
to obtain the catalogue of radio sources. To estimate the efi'ect of the resolution bias on our observations we have 
modelled the effective angular-size distribution of the sources in our sample and we have used this distribution to 
simulate a sample of radio sources. Finally we present the radio count distribution down to 0.08 mjy derived from 
the catalogue. Our counts are in good agreement with the best fit derived from earlier surveys, and are about 50% 
higher than the counts in the HDF. The radio count distribution clearly shows, with extremely good statistics, 
the change in the slope for the sub-mjy radio sources. 
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1. Introduction 

It is well established that the 1.4 GHz source counts 
at sub-mJy levels reveal the presence of a popula- 
tion of faint radio sources far in excess with respect 
to those expected from the high luminosity radio 
galaxies and quasars which dominate at higher fluxes 
(jWindhorst et al. 1985irCondon 1989: ^ Hopkins et al. 1998| 
Ciiiegi et ai. 1999| rRIcEards 20001 IPrimdoni et al. 2001al 
nnippion^^L^1999b| ) . Early spectroscopic stud- 
ies, limited to relatively bright optical counterparts 
(B<21. 5-22.0; IBenn et al.T993|l . suggested that most 
of these sub-mJy radio sources were starburst galaxies. 
However it has been shown that the predominance of 



starburst galaxies is dependent on the magnitude limit 
of the spectroscopic follow up ( |Gruppioni et al. 1999aj 
IPrandoni et al. 2001b|l . While at bright magnitude 
(B ^ 22) most of the optical counterparts are indeed star- 
burst galaxies, at fainter magnitudes (B>22.5) most of 
the optical counterparts appear to be early type galaxies. 
This mixture of at least two different populations is 
consistent with what is being found at even fainter radio 
fluxes, in the /iJy regime, where high-z early type galax- 
ies, intermediate-z post starburst galaxies, and lower-z 
emission line galaxies are found in approximately similar 
proportions IjHammer et al. 19951 IWindhorst et al. 19951 
IRichards et al. 1998|l . 
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Fig. 1. South-eastern quadrant of the VLA-VDF radio image. Contours are in units of signal to noise ratio, first 
contour is 5cr. 



In order to fully investigate the nature and evolution 
of the sub-mJy population it is absolutely necessary to 
couple deep radio and optical (both imaging and spectro- 
scopic) observations over a reasonably large area of the 
sky. 

The VIRMOS VLT Deep Survey (VVDS, 
ILe Fevre et al. 2002|l will produce spectroscopic red- 
shifts for about 1.5 x 10^ galaxies in an area of ~ 16 
deg^ selected from an unbiased photometric sample of 
more than 1 million galaxies. We have selected a 1 deg^ 
field from the VVDS, centered at RA(J2000)=02:26:00 
DEC(J2000)=-04:30:00, for deep VLA radio observa- 
tions at 1.4 GHz (hereafter the VLA-VIRMOS Deep 
Field, VLA-VDF). This field is ideal for a radio survey 
as UBVRI photometry, complemented by K band data 
on a smaller region, is already available to Iab — 25 
IjLe Fevre et al. 2 ' QOI ) , and spectroscopy is being obtained 
to Iab = 24 - 24.5 with the VIMOS spectrograph at the 
VLT ULe Fevre et al. 2nni|l . 

In this paper we present the VLA radio observations 
at 1.4 GHz of the VLA-VIRMOS Deep Field, discuss the 
methods used to derive the catologue of about 1000 radio 
sources (down to a limit of 80 /iJy) and derive the radio 
source counts. The identification of the radio sources using 
the multi-band photometry will be discussed in a follow- 
ing paper ( |Ciliegi et al. 200'2| ). In Sect. 2 the observations 
and data reduction are presented. Sect. 3 contains a de- 
tailed description of the analysis carried out on the radio 
mosaic in order to quantify the effects of clean bias and 



bandwidth smearing on our observations. The procedure 
adopted to obtain a complete catalogue of radio sources 
from the radio mosaic is presented in Sect. 4. Finally, in 
Sect. 5 we derive the radio counts corrected for the reso- 
lution bias. Conclusions are given in Sect. 6. 

2. The VLA Observations 

The observations were obtained at the Very Large Array 
(VLA) in B-configuration for a total time of 56 hours over 
9 days from November 1999 to January 2000. This config- 
uration was adopted as the best choice in order to obtain 
a deep survey of a relatively large (1 square degree) area 
of the sky with an acceptable resolution. 

At 1.4 GHz the VLA antennas have a primary beam 
with a FWHM of 31 arcmin. In order to image with uni- 
form sensitivity a 1 square degree field it is necessary to 
make multiple pointings displaced by about 31/-\/2 ~ 22 
arcmin HCondon et al. 1998llBecker et al. 1995| . We chose 
to cover the surveyed area with a square grid of 9 point- 
ings, separated by 23 arcmin in right ascension and decli- 
nation. Such a geometry allows to reach theoretical noise 
variations smaller than 10% over 95% of the 1 deg^ field. 
Each of the pointing centers was observed for a total 
of about 6 hours, including the observations of the cal- 
ibrators. Every 20 minutes we interleaved the scans on 
the nine pointings with a short observation of the source 
J0241— 082 to provide amplitude, phase, and bandpass cal- 
ibration. 
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Fig. 2. South-western quadrant of the VLA-VDF radio image. Contours are in units of signal to noise ratio, first 
contour is 5f7. 



The observations were carried out in bandwidth syn- 
thesis mode to avoid substantial chromatic aberration 
(bandwidth smearing). In this way it is also possible to 
reduce the effects of narrow-band interferences since only 
the channel affected by the interferences, instead of the 
whole bandwidth, can be removed from the data. The 
data were collected in spectral line mode using two in- 
termediate frequency (IF) bands centered at 1364.9 MHz 
and 1435.1 MHz. Each IF was divided in 7 channels each 3 
MHz wide. Due to limitations in the VLA correlator only 
circular polarization modes were recorded. 

2.1. Calibration & Editing 

The data were reduced and analyzed using the pack- 
age AIPS developed by the National Radio Astronomy 
Observatory. The amplitude calibration was derived from 
daily observations of 3C 48 assuming flux densities of 16.51 
Jy at 1365 MHz and 15.87 Jy at 1435 MHz. The AIPS 
tasks UVLIN and CLIP were used to flag bad visibility 
data resulting from radio frequency interferences, receiver 
problems or correlator failures. The task UVLIN was run 
after the amplitude calibration with a threshold of 1 Jy. 
After UVLIN, the task CLIP was run on a channel-by- 
channel basis with a threshold determined by the task 
UVPRM (8 times the rms of the data). The task UVFIX 
was run to improve the astrometric solution and then the 
data sets for each pointing were combined and inspected 
for further editing. 



2.2. Imaging, self calibration and mosa icing 

Self calibration and imaging of wide field deep observa- 
tions is a time consuming task. For each pointing we im- 
aged a 2048 X 2048 pixels area (51 x 51 arcminutes, 1 pixel 
corresponds to 1.5 arcsec) along with a number of smaller 
images (usually 32 x 32 pixels) centered on off-axis sources 
that can produce confusing grating rings in the imaged 
area. The possibly confusing sources have been identified 
with the RUN file generator applet available at the NVSS 
home page, selecting all the sources with peak flux den- 
sity greater than 1 mJy (at the NVSS resolution) within 
a radius of 60 arcmin of each pointing and not included 
in the main field area. 

To avoid distorsions due to the use of two dimen- 
sional FFT to approximate the curved celestial sphere, the 
2048 X 2048 pixels area of each pointing was not decon- 
volved as a single image but was split in a number of sub- 
images (e.g. |Perley 1999| ). At the end of the self-calibration 
deconvolution iteration scheme, the sub- images were com- 
bined together using the AIPS task PASTE to produce the 
final 2048 x 2048 pixels image of each single pointing. The 
final images have been restored with a 6 x 6 arcseconds 
FWHM gaussian beam. We self-calibrated and cleaned the 
different pointings in a way as homogeneous as possible in 
order to minimize differences in the sensitivity. Clearly, the 
presence or absence of relatively strong (~ 10 mJy/beam) 
sources in some fields and the fluctuations in the noise pro- 
duced slightly different noise figures, with 1 a rms noise 
ranging from 14.8 to 17.9 /iJy/beam, in the nine pointings. 
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Fig. 3. North-eastern quadrant of the VLA-VDF radio image. Contours are in units of signal to noise ratio, first 
contour is 5f7. 



Finally, the 9 pointings have been combined together using 
the task HGEOM and LTESS obtaining a linear combi- 
nation weighted by the square of the beam response. The 
average noise over the full 1 square degree field in the mo- 
saic map is f 7.5 /iJy. In Fig. lfl4l we show the final image of 
the 1 square degree VLA-VIRMOS deep field split in four 
quadrants for a clearer representation. The noise over the 
1 square degree field is homogeneous (see also Sect. 4.1) 
and the few regions devoid of sources visible in Fig. 1 1141 
(the most notable of which is the area around right ascen- 
sion 02:25:40 and declination -04:52:00) are real and not 
artifacts produced by a much higher local noise. 

3. Analysis 

3.1. Clean bias 

Clean bias is a recently recognized possible source of er- 
ror in the flux density estimate derived from interferom- 
eter snapshot observations. The FIRST and NVSS VLA 
surveys, for example, are affected by this bias and the 
flux densities derived from the publicly available images 
have to be corrected a-posteriori. The effect of the clean 
bias is that flux densities of real sources are sistemati- 
cally underestimated because the CLEANing algorithm 
subtracts flux from real sources and redistributes it on 
top of noise peaks or sidelobes. This effect is dependent 
on the image noise and synthesized beam sidelobe lev- 
els and mostly independent from the source flux density 
^Recker et a.l. 1995..Condon et al. IflflSjl. 



The VLA-VDF observations are very long compared 
to the snapshots of the VLA surveys, and consequently 
the synthesized beam has much lower sidelobes. For in- 
stance, the NVSS observations have a synthesized beam 
with sidelobes reaching about 25% of the main lobe, while 
the sidelobes in the synthesized beam of the VLA-VDF ob- 
servations reach only 1.6% of the main lobe. A good rule 
of thumb to avoid clean bias is to minimize the clean area 
and not to clean the image down to the theoretical noise. 
The first prescription is hard to follow as our goal is to 
image a 1 square degree field. On the other hand, we used 
a rather conservative approach halting the clean process 
when the clean residuals were between 2 and 5 times the 
theoretical r.m.s. noise. While we can expect that the clean 
bias in our observations is much lower than that affecting 
snapshot observations, we can not rule it out completely 
a-priori. In order to assess the impact of the clean bias 
on the VLA-VDF observations we have used two different 
methods. 

The first method is a step-by-step clean. We can expect 
that beyond a threshold value in the number of iterations, 
the clean bias begins to be important and the flux density 
of the components in these maps becomes sistematically 
lower than that of the corresponding components in the 
images with less iterations. About 4 x 10'* clean iterations 
are the minimum number necessary to effectively approach 
the expected noise for the 2048 x 2048 pixel images. We 
have then chosen one of the nine fields and produced differ- 
ent images with an increasing number of clean iterations. 
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Fig. 4. North-western quadrant of the VLA-VDF radio image. Contours are in units of signal to noise ratio, first 
contour is 5cr. 



Table 1. Clean bias simulation 



Niter (1 


-Niter (2) 


median AS 


40000 


- 50000 


0.5 ^Jy/beam 


50000 


- 75000 


1.0 ^Jy/beam 


75000 


- 500000 


14.4 ^Jy/beam 



We have produced images with 4 x 10^, 5 x 10^, 7.5 x 10"* 
and 5 x 10'"' iterations and on each of them we have identi- 
fied radio components down to the 3cr limit. For each com- 
ponent we have extracted the peak flux and computed the 
difference of fluxes obtained between maps with different 
iteration limits. The median of these differences is listed 
in Table ^ As can be seen from Table increasing the 
number of iterations to 5 x 10^ or even 7.5 x 10"* produces 
a systematic decrease in the peak flux densities of about 
1.0 ^Jy/beam or less. We have to push the clean to a far 
higher number of iterations to start seeing a significant ef- 
fect due to the clean bias. Using 5 x 10^ iterations we find 
a decrease in the peak flux density of the sources with a 
median of 14.4 /iJy/beam. We have then chosen a limit of 
about 5 X lO'* clean iterations for each of the nine fields. 
This means that any clean bias affecting our brightness 
measurements should be less than 1 /xJy/beam and then 
practically negligible. 

The second method used to verify the absence of a 
significant effect produced by the clean bias was to insert 
artificial sources with known flux in the uv-data set of a 
chosen field. In particular, we have modified the uv-d&ta 



set of a randomly chosen field adding 25 point sources with 
0.5 mJy flux. We have then cleaned the field to the same 
depth used for the original one and compared the fluxes 
of the artificial sources on the map with their true fluxes. 
We have repeated this operation four times obtaining a set 
of 100 artificial sources at different positions. The mean 
of the peak flux density distribution of these 100 sources 
is 0.503 ± 0.003 mJy/beam. We can conclude that both 
methods confirm that the flux density derived from our 
images are not affected by the clean bias. 

3.2. Bandwidth Smearing 

Imaging sources at large distances from the phase center 
can result in radial smearing reducing the peak flux den- 
sity of a source while conserving its integrated flux density. 
This effect is known as bandwith smearing (or chromatic 
aberration) and affects all the synthesis observations made 
with a finite bandwidth. The image smearing is propor- 
tional to the bandwidth and to the distance of the source 
from the phase center. In order to image a 1 square degree 
field we have to minimize the effect of bandwith smear- 
ing and for this reason we observed in spectral line mode. 
Nonetheless, some amount of smearing can still be present 
in our images. To check this effect on our observations we 
have observed the radio source 3C84 at different position 
offsets (0, 5, 10, 20, 30 arcmin) in two orthogonal direc- 
tions. The images of 3C84 at different offsets have been 
fitted with a two dimensional gaussian to derive the peak 



6 



M. Bondi et al.: The VLA-VIRMOS Deep Field 




10 20 
Off — axis (arcmin) 

Fig. 5. Ratio between Sp and 5t for 3C 84 as a function of 
off-axis angular distance as a test of bandwidth smearing 
for the VLA-VDF observations. The vertical line indicates 
the distance between adjacent pointing centers. 

and integrated flux density. The mean ratio between these 
two quantities is shown in Fig. At the distance between 
the pointing centers (23 arcmin) Sp/St ~ 0.96. Since the 
nine radio fields are combined together weighted by the 
square of the beam response, we can conclude that the 
effect of the bandwidth smearing is negligible. 

4. From the radio mosaic to the catalogue 

4.1. The noise map 

Having verified that clean bias and radial smearing do 
not significantly affect the determination of the flux den- 
sities in the VLA-VDF, we started to work on the ex- 
traction of a catalogue from the mosaic image. In or- 
der to select a sample above a given threshold, defined 
in terms of local signal to noise ratio, we performed a 
detailed analysis of the spatial root mean square (rms) 
noise distribution over the entire mosaic image using the 
software package SExtractor (IBertin fc Arnouts 1996|l . To 
construct the noise map, SExtractor makes a first pass 
through the pixel data, computing an estimator for the 
local background in each mesh of a grid that covers the 
whole frame fsee lBertin &: Arnouts 1996l for more details). 
The choice of the mesh size is very important. When it is 
too small, the background estimation is affected by the 
presence of real sources. When it is too large, it cannot 
reproduce the small scale variations of the background. 
For our radio mosaic we adopted a mesh size of 20 pix- 
els, corresponding to 30 arcsec. A grey scale of the noise 
map obtained with SExtractor is shown in Fig. (HJ while 
Fig. [3 shows the histogram of its pixel values. Due to the 
uniform noise over the whole field the pattern of the 9 



Right A3ceri5ir;n (ZOQD) 



Fig. 6. A grey scale of the noise map (I x 1 deg^) obtained 
with SExtractor. Darker regions mean higher noise. 
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Fig. 7. Distribution of the pixel values of the noise map 
obtained with SExtractor 



pointings used for our observations can be barely seen in 
Fig. ini The areas of higher noise (black pixels) are due 
to the presence of relatively strong radio sources (10-20 
mjy/beam) in the map. The pixel values distribution has 
a peak at 16 /iJy/beam, well in agreement with the noise 
values found in Sect. 3.2. 

Since SExtractor was developed for the analysis of op- 
tical data and this is the first attempt in using it on a radio 
map, we tested the reliability of its output by construct- 
ing a noise map with a completely independent software 
written in IDL language. Briefly, starting from the resid- 
ual map obtained from the AIPS task SAD where all the 
sources with peak flux greater than 60 /iJy (^ 3.5cr) have 
been subtracted, we have first removed the most anoma- 
lous residual pixels (including extended sources not found 
or rejected by SAD) substituting these values with the 
average rms obtained from the residual map. Then, we 
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applied a local sigma clipping, substituting all the pixel 
values greater than 3 times the local sigma with a random 
value extracted from a Gaussian with mean and sigma 
equal to the local values. Finally, the noise map has been 
obtained substituting each pixel value with the standard 
deviation values calculated in a local box around each 
pixel (we used a box of 20x20 pixels). In Fig.|51we show 
the distribution of the ratio (pixel by pixel) between the 
noise map obtained with SExtractor and that obtained 
with our IDL code. A detailed analysis of the two maps 
shows that they are in very good agreement with each 
other, with differences of pixel values which are smaller 
than 25% over about 96% per cent of the area. 

On the basis of this comparison we concluded that the 
noise map obtained with SExtractor is indeed reliable and 
we used it for the extraction of a catalogue. 

4.2. The source detections 

The area from which we have extracted the complete 
sample of sources is 1x1 deg^ centered at RA=02:26:00 
DEC=-04:30:00 (J2000). Within this region we extracted 
all the radio components with a peak flux Sp >60 /iJy 
(~ 3.5ct) using the AIPS task SAD (Search And Destroy), 
which attempts to find all the components whose peaks 
are brighter than a given level. For each selected compo- 
nent, the peak and total fluxes, the position and the size 
are estimated using a Gaussian fit. However, for faint com- 
ponents the Gaussian fit may be unreliable and a better 
estimate of the peak flux Sp (used for the selection) and 
of the component position is obtained with a simple in- 
terpolation of the peak values around the fitted position. 
Therefore, starting from the SAD positions, we derived the 
peak flux Sp and the position of all the components using 
a second-degree interpolation with task MAXFIT. Only 
the components for which the ratio between the MAXFIT 
peak flux density and the local noise (derived from the 
noise map described in the previous section) was greater 



or equal to 5 have been included in the sample. A total of 
1084 components have been selected with this procedure. 

Some of these components clearly belong to a single 
radio source (e.g. the lobes of a few FR II radio sources, 
or components in very extended sources), but for most of 
them it is necessary to derive a criterion as general as pos- 
sible to discriminate between different components of the 
same radio source or truly different radio sources. For this 
reason, the components with distance less than 18 arcsec- 
onds (3 times the beam size) have been selected as possible 
doubles and have been visually checked one-by one using 
also preliminary deep optical images. Based on the com- 
parison of the radio and optical fields, we have assumed 
that when two components have a distance smaller than 18 
arcsec, a flux ratio smaller than 3, and both components 
have peak brightness larger than 0.4 mJy/beam they be- 
long to the same radio source, otherwise are considered as 
separate radio sources. With this choice, considering the 
number of sources with flux greater than 0.4 mJy we can 
expect 3 spurious couples of radio components in the 1 
deg^ field, compared with the 40 observed. On the other 
hand, within 18 arcsec we expect 50 and find 51 couples 
of components with both fluxes less than 0.4 mJy. 

The final catalogue lists 1054 radio sources, 19 of which 
are considered as multiple, i.e. fitted with at least two 
separate components, and it will be available on the web at 
http://virmos.bo.astro.it/radio/catalogue.html A sample 
page of the catalogue is shown in Table HJ 

For each source we list the source name, position in 
RA and DEC with errors, peak flux and total flux density 
with errors, major and minor axis and position angle. For 
the unresolved sources the total flux density is equal to the 
peak brightness and the angular size is undetermined. For 
each of the 19 sources fitted with multiple components we 
list in the catalogue an entry for each of the components, 
identified with a trailing letter (A, B, C, . . . ) in the source 
name, and an entry for the whole source, identified with a 
trailing T in the source name. In these cases the total fiux 
was calculated using the task TVSTAT, which allows the 
integration of the map values over irregular areas, and the 
sizes are the largest angular sizes. The peak fiux density 
distribution of the 1054 radio sources is shown in Fig. O 

4.3. Resolved and Unresolved Sources 

Since the ratio between total and peak fluxes is a direct 
measure of the extension of a radio source, we used it 
to discriminate between resolved or extended sources (i.e. 
larger than the beam) and unresolved sources. 

In Fig. ^1 we plot the ratio between the total St and 
the peak Sp flux density as a function of the peak flux 
density for all the radio sources in the catalogue. To select 
the resolved sources, we have determined the lower enve- 
lope of the fiux ratio distribution of Fig- ITUI and. assuming 
that values of Sx/Sp smaller than 1 are purely due to sta- 
tistical errors, we have mirrored it above the St/ Sp = l 
value (upper envelope in Fig. llO|l . We have considered ex- 
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24 
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Fig. 9. Peak flux density distribution for all the 1054 
VIRMOS radio sources 



Fig. 11. Visibility area of the VIRMOS radio survey. It 
represents the area over which a source with a peak flux 
density Sp can be detected. 
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Fig. 10. Ratio of the total flux St to the peak flux Sp 
as a function of Sp. The solid lines show the upper and 
lower envelopes of the flux ratio distribution containing 
the sources considered unresolved (see text). Open simbols 
show the sources considered extended. 



tended the 254 sources laying above the upper envelope, 
that can be characterized by the equation 



0.95-(i/Sp) 



(1) 



4.4. Errors in the source parameters 



The formal relative errors determined by a Gaussian fit 
are generally smaller than the true uncertainties of the 
source parameters. Gaussian random noise often domi- 
nates the errors in the data l|Condon 1997|l . Thus, we used 
the Condon (1997) error propagation equations to esti- 
mate the true errors on fluxes and positions: 



SI 



(2) 



where Sp and St are the peak and the total fluxes, 



and p is the signal-to-noise ratio, given by 



4^2 



map 



(3) 



are the fitted FWHMs of the major 
map is the noise variance of the image 



where 9m and 
and minor axes, a- 
and is the FWHM of the Gaussian correlation length 
of the image noise (~FWIIM of the synthesized beam). 

The exponents are om = am = 3/2. 

The projection of the major and minor axis errors onto 
the right ascension and declination axes produces the total 
rms position errors given bv lCondon et al. 19981 



= £i + < sin^(PA)+a^„ cos' (PA) 
= + < cos\PA) + sin' [PA) 



(4) 
(5) 



where PA is the position angle of the major axis, 
{Sa,s&) are the 'calibration" errors, while CTxo and ay^ are 
el^/{Aln2)p' and 9l,/{4ln2)p' respectively. 

Calibration terms are in general estimated from com- 
parison with external data of better accuracy than the one 
tested, using sources strong enough that the noise terms 
in equation ^ and |S1 are much smaller than the calibra- 
tion terms. Unfortunately there are no such data available 
in the region covered by the VLA-VDF survey. In fact, 
the only other radio data available in this region are the 
NVSS radio data, with a synthesized beam of 45", about 
a factor 8 greater than the synthesized beam of the VLA- 
VDF survey. To estimate the calibration terms and 
es we have selected all the point sources with Sp > 0.3 
mJy/beam from the final mosaic (105 objects) and com- 
pared their positions with those found on the single im- 
ages. The mean values and standard deviations found from 
this comparison are < ARA>= 0.007 ± 0.076 arcsec and 
< ADEC>= -0.039 ±0.115 arcsec. These values are con- 
sistent with no systematic offset in right ascension and 
declination and a standard deviation of about 0.1 arcsec. 
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Fig. 12. Integral angular size distribution (solid line) for 
the sources in the VLA-VDF survey with 0.4<S<1.0 mJy. 
The dashed line shows the best fit obtained for this distri- 
bution (h(6l)=l/(1.6^) for < 4" and h(0)= {9-^-^-0.01) 
for 9 > 4"), while the dot-dashed line shows the integral 
angular distribution reported by Windhorst et al. (1990) 
and assumed by several authors as the true angular size 
distribution of radio sources. 



In calculating the errors affecting the radio position of the 
sources in the catalogue we have assumed Eq, ~ ~ 0.1 
arcsec. 



5. Survey Completeness and Source Counts 

The visibility area of the VIRMOS radio survey as a func- 
tion of the peak flux density Sp is shown in Fig. ^] As 
expected, the visibility area increases very rapidly between 
0.05 and 0.09 mJy and becomes equal to 1 degree at Sp > 
0.093 mJy. This is a consequence of the observing strat- 
egy that has assured a very uniform noise over almost the 
entire 1 square degree field used for the extraction of the 
catalogue. Since the completeness of the radio sample is 
defined in terms of the peak flux, while the source counts 
will be derived as a function of the total integrated flux, 
corrections must be applied to the observed numbers of 
radio sources in order to take into account all possible ob- 
servational biases. The most important of such biases is 
probably the resolution bias which leads to missing faint 
(i.e extended) sources at fluxes close to the limit. In fact, 
such sources, with peak flux densities below the survey 
limit, but total integrated fluxes above this limit, would 
not appear in the catalogue. The correction due to this 
bias is a function of the intrinsic angular size distribu- 
tion of the sources and of the beam of the observations. 
To estimate the correction factor to be applied to the ob- 
served data, in the next section we will model the effective 
angular-size distribution of the sources in our radio sample 
and then we will use this distribution to simulate a sample 
of radio sources that we will analyse with the same recipe 
used for the real sources described in Sect. 4. 



5.1. Angular Size Distribution 

Previous high-resolution studies of the faint radio pop- 
ulation suggested that the median angular size {9med) 
for sub-mjy radio sources is approximately 2" and al- 
most independent of flux density between 0.08-1.0 mJy 
IjWiTidhorst et al. 1 998llTW^Jont et al. 1991ll()ort 198811. 
In order to derive an unbiased distribution of angular sizes 
from our sample, we have to use sources in a range of to- 
tal fluxes in which the resolution bias is likely not to have 
modifled, in the catalogue, the intrinsic angular size dis- 
tribution. For this purpose we used all the sources with 
0.4 < S < 1 mJy. Given the relation between angu- 
lar size 9 and the ratio between total and peak fluxes, a 
source with St ~ 0.6 mJy would have a peak flux greater 
than our detection limit even for a relatively large angu- 
lar size {6 ~ 15 arcsec). Forty-eight of the 111 sources 
(43%) in this flux range are resolved, and we measured 
for them a value of 9. The other 63 sources are unre- 
solved, i.e. they are below the solid line representing the 
upper envelope in Fig 1101 For them we could derive only 
an upper limit to their intrinsic size 9] these upper lim- 
its range from ~ 1.4 to ~ 2.3 arcsec. These 9 values 
were then analyzed with the survival analysis techniques 
of iFeigelson & Nelson (1985)' using the statistical pack- 
age ASURV (Isobc ct al. 1986). This technique uses all 
the upper limits (which are slightly more than 50% in our 
data set of angular sizes) in reconstructing the intrinsic 
distribution. The resulting estimate for the median value 
of 9 is ^mcd — 1.8 ± 0.2 arcsec, somewhat lower than, but 
consistent with the value ^med — 2.6 ± 0.4 obtained by 
Richards (2000) for sources in the same flux range from 
his deep survey in the Hubble Deep Field region. Our 
value of ^med is also in good agreement with the rela- 
tion 9^ed{") = 2.0S"-^, found by Windhors fet al 19901 
where S is the total flux in mJy. In Fig. 1121 we report 
the integral angular size distribution for our sources with 
0.4<S<1.0 mJy (solid line) which, because of the pres- 
ence of the upper hmits, is determined only for 9 > 1".8. 
The dashed line shows an analytical flt to this distribu- 
tion (h(6l)=l/(1.6^) for 9 < 4" and h(6')= {9-^'^ - 0.01) 
for 9 > 4"). For comparison the dot-dashed line shows the 
integral angular distribution h(0)=exp[— In 2(0/0med)°'^^] 
with 9yncd = 1.8 suggested bv Windhorst et al. 19901 

As clearly shown in Fig.^J the Windhorst et al. rela- 
tion is not a good representation of our measured distribu- 
tion of angular sizes, because it predicts a substantial tail 
of sources with large angular sizes, which is not present in 
our data. For example, the fraction of sources with 9 > A 
arcsec predicted by this relation (~ 32%) is about twice 
as high as that measured from our data, as well as that 
which can be derived for sources in the same flux range 
from the radio data in the HDF region fRich ards 20001 see 
his Fig. 4). The use of the Windhorst relation, because of 
its high fraction of sources with large angular sizes, would 
lead to an overstimate of the correction factors due the 
resolution bias. 
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S 


< S > 


Ns 


dN/dS 




C 


iV(> S) 
deg-^ 


(mJy) 


(mJy) 




sr-i Jyi-5 






0.08 


- 0.12 


0.10 


346 


3.07x10^" 


2.91± 0.16 


1 


.29 


1200 


0.12 


- 0.18 


0.15 


205 


1.12x10^° 


2.94± 0.21 


1 


.25 


718 


0.18 


- 0.27 


0.22 


167 


6.09x10" 


4.40± 0.34 


1 


.00 


462 


0.27 


- 0.41 


0.33 


99 


2.41x10" 


4.79± 0.48 


1, 


.00 


295 


0.41 


- 0.61 


0.50 


58 


9.40x10* 


5.15± 0.68 


1 


.00 


196 


0.61 


- 0.91 


0.74 


42 


4.54x10* 


6.85± 1.06 


1 


.00 


138 


0.91 


- 1.37 


1.12 


25 


1.80x10* 


7.50± 1.50 


1 


.00 


96 


1.37 


- 2.05 


1.67 


18 


8 65x10^ 


9.91± 2.34 


1 


.00 


71 


2.05 


- 3.08 


2.51 


17 


5.44x10^ 


17.20± 4.17 


1 


.00 


53 


3.08 


- 4.61 


3.77 


12 


2.56x10^ 


22.31± 6.44 


1 


.00 


36 


4.61 


- 6.92 


5.65 


6 


8.54x10^ 


20.49± 8.37 


1 


.00 


24 


6.92 


- 10.38 


8.48 


7 


6.64x10^ 


43.92±16.60 


1 


.00 


18 


10.38 


- 15.57 


12.71 


4 


2.53x10*^ 


46.10±23.05 


1 


.00 


11 


15.57 


- 119.23 


42.91 


7 


2.24x10^ 


85.32±32.24 


1 


.00 
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5.2. Completeness 



In order to estimate the combined effects of noise, source 
extraction and flux determination techniques and reso- 
lution bias on the completeness of our sample, we con- 
structed simulated samples of radio sources down to a flux 
level of 0.04 mJy, i.e. a factor of two lower than the min- 
imum flux we used to derive the source counts (see Sect. 
5.3). This allows to take into account also those sources 
with an intrinsic flux below the detection limit which, be- 
cause of positive noise fluctuations, might have a mea- 
sured flux above the limit. The simulated samples have 
been extracted from source counts with the integral size 
distribution derived in the previous section and described 
by a broken power law consistent with that observed. 



Following these recipes, we simulated 9 samples, each 
of them with a number of sources above the detection limit 
similar to that observed in the real data (i.e. ~ 1000). 
All the sources, including those below the detection limit, 
were randomly injected in the CLEANed sky images of 
the field and were recovered from the image and their 
fluxes were measured using the same procedures adopted 
for the real sources (see Sect. 4). The detected simulated 
sources were then binned in flux intervals. Finally, from 
the comparison between the number of simulated sources 
detected in each bin and the number of sources in the in- 
put simulated sample in the same flux bin we computed 
the correction factor C to be applied to our observed 
source counts. In Table |31 we report the average correc- 
tion factor C for each flux density bin. As expected, the 
resolution bias significantly affects the first two flux den- 
sity bins in the source counts. Our simulations tell us we 
are missing about 29% and 25% of sources in the first two 
flux density bins respectively. In the bins at higher flux the 
results of these simulations are consistent with no need for 
a correction. We therefore set C = 1.0 for all these bins. 



5.3. Source counts 

In order to reduce possible problems near the flux limit 
and to avoid somewhat uncertain corrections in the steep 
part of the visibility area we constructed the radio sources 
counts considering only the 1013 sources with a flux den- 
sity greater than 0.08 mJy {i.e. we excluded the 41 sources 
with S<0.08 mJy). The 1.4 GHz source counts of our sur- 
vey are summarised in Table 13 For each flux density bin, 
the average flux density, the observed number of sources, 
the differential source density (in sr"-'^ Jy~^)7 the nor- 
malised differential counts nS'^'^ (in sr~^ Jy^'^) with es- 
timated Poissonian errors (as n^/^S*^'^). In the last two 
columns we report the correction factor C to be applied 
to our source counts to correct for incompleteness and the 
corrected integral counts (in deg^^). 

The normalised 1.4 GHz counts (column V) multiplied 
by the correction factor C are plotted in Fig.ll3lwhere. for 
comparison, the differential source counts obtained with 
other 1.4 GHz radio surveys are also plotted. 

Our counts are in good agreement, over the entire flux 
range sampled by our data (0.08 - 10 mJy), with the best 
fit to earlier surveys ( |Katgert et al. 1 988). It is interesting 
to note that, because of the relatively good statistics of our 
data points over about two orders of magnitude in flux, 
our data clearly show the change in slope of the differen- 
tial counts, occurring below 1 mJy. Fitting the VLA-VDF 
differential and integral counts with two power laws we 
obtain, for S < 0.6 mJy: 

^ = (57.54 ±1.07) X 5-2.28±o.o4 
as 

and at higher fluxes, S > 0.6 mJy: 

^ = (75.86 ± 1.08) x5~i ™±0 "5 
as 

where dn/ds is in deg^^mJy"^ e 5 in mJy. At faint 
fluxes {S < 0.4 mJy) the most relevant comparison is 
with the results obtained by Richards (2000) in his deep 
VLA survey of a region of about 40 x 40 arcmin around 
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Fig. 13. The normalised differential source counts. The abscissa gives the flux density (mJy) and the ordinate gives the 
differential number of sources normalised by S'^'^ (sr~^ Jy"*"^'^). With different symbols are reported source counts at 
1.4 GHz from several surveys: the smaller open circles represent the counts from the FIRST survey l| White et a. 1997|l . 
the open squares are the counts from the ATESP survey ('Prando ni et al. 2001a|l . the larger open circles are the counts 
obtained by the survey ELAIS ( |Gruppioni et al. 1999b[Ciliegi et al. 1999| ), the stars are the counts from the HDF 
l|R,ichards 2000|l . the asterisks are the counts from the Phoenix Deep Survey ( [Hopkins et al. 1998| ), the filled triangles 
are the counts obtained with the VIRMOS survey discussed in this paper. The solid lines is a best fit to a compilation 
of earlier deep surveys made by .Katgert et al. 19'88| 



the Hubble Deep Field. In this flux range the statistics 
in our data (i.e. number of radio sources) is about four 
times higher than that in the HDF data and our derived 
counts are about 50% higher than those of Richards (see 
also [Hopkins et al. 2002| for similar results). This is con- 
sistent with the fact, already noted by Richards, that the 
counts in the HDF region appear to be sistematically lower 
than those of other fields above 0.1 mJy. This difference 
is probably due to a combination of real field-to field vari- 
ations and some possible incompleteness of somewhat ex- 
tended sources in the high resolution HDF data (2" beam). 
Indeed, lower resolution, deep Wcsterbork observations of 
a sub-area in the Hubble Deep Field region have detected 
a non-negligible number of sources which do not appear 
in Richards' catalogue ([Garrett et al. 2000(1 . 



6. Conclusions 

Using the VLA at 1.4 GHz we have observed a 1 deg^ 
field centered on the VIRMOS Deep Field (RA=02:26:00 
DEG=— 04:30:00), imaging the whole area with uniform 
sensitivity and a resolution of 6 arcsec. We have investi- 
gated the effects of clean bias and bandwidth smearing 
on our observations confirming that the observing strat- 
egy and the data reduction procedure allow us to con- 
sider these effects negligible. A complete catalogue of 1054 
radio sources down to a local 5(7 limit (~ 80 ^Jy) has 
been compiled. In order to assess the effects of random 
noise, source extraction technique and resolution bias on 
the completeness of our sample we have first derived the 
effective angular-size distribution using the sources in our 
sample in the range 0.4-1.0 mJy. Then we have gener- 
ated a large sample of simulated sources using the derived 
angular-size distribution and extracted from source counts 
described by a broken power law consistent with that ob- 
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served. These simulations allowed us to statistically cor- 
rect ourc counts in the faintest flux bins. The final counts 
are in good agreement with the best fit to earlier surveys 
UKatgert et al. 1988| ). In particular, our data clearly show 
a significant change in slope of the differential counts oc- 
curring below 1.0 mJy. The best fit slope in the range 0.08- 
0.6 mJy (-2.28 ±0.04) is close to the Euclidean value. At 
faint fluxes {S < 0.4 mJy), where we have a high statis- 
tics, our derived counts are about 50% higher than those 
of IRichards 20T)f)l in the HDF region. This is consistent 
with the fact, already noted by Richards, that the counts 
in the HDF region appear to be sistematically lower than 
those of other fields above 0.1 mJy 

The same region of the sky has been target of deep 
(-^AB = 25), multicolor (UBVRIK) photometry observa- 
tions. The photometric identification of most of the radio 
sources in the catalogue and planned spectroscopic obser- 
vations during the VIRMOS Deep Field Survey will pro- 
vide a unique opportunity to study the nature and prop- 
erties of the fiJy source population. 
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